arXiv: 1 402 . 1 5 3 8 v 1 [astro-ph. SR] 7Feb2014 


High-Resolution Submillimeter and Near-Infrared Studies of the 

Transition Disk around Sz 91 

Takashi Tsukagoshi 1 , Munetake Momose 1 , Jun Hashimoto 5 , Tomoyuki Kudo 2 , Sean 
Andrews 4 , Masao Saito 2 , Yoshimi Kitamura 3 , Nagayoshi Ohashi 2 , David Wilner 4 , Ryohei 
Kawabe 2 , Lyu Abe 7 , Eiji Akiyama 2 , Wolfgang Brandner 8 , Timothy D. Brandt 9 , Joseph 
Carson 10 , Tliayne Currie 28 , Sebastian E. Egner 11 , Miwa Goto 12 , Carol Grady 13 , Olivier 
Guyon 11 , Yutaka Hayano 11 , Masahiko Hayashi 2 , Saeko Hayashi 11 , Thomas Henning 8 , Klaus 
W. Hodapp 14 , Miki Ishii 2 , Masanori lye 2 , Markus Janson 15 , Ryo Kandori 2 , Gillian R. 

Knapp 15 , Nobuhiko Kusakabe 2 , Masayuki Kuzuhara 2,16 , Jungmi Kwon 27 , Mike 
McElwain 17 , Taro Matsuo 18 , Satoshi Mayama 27 , Shoken Miyama 19 , Jun-ichi Morino 2 , 
Amaya Moro-Martin 20 , Tetsuro Nishimura 11 , Tae-Soo Pyo 11 , Eugene Serabyn 21 , Takuya 
Suenaga 27 , Hiroshi Suto 2 , Ryuji Suzuki 2 , Yasuhiro Takahashi 6 , Hideki Takami 2 , Michihiro 
Takami 22 , Naruhisa Takato 11 , Hiroshi Terada 11 , Christian Thalmann 23 , Daigo Tornono 11 , 
Edwin L. Turner 9 , Tomonori Usuda 11 , Makoto Watanabe 24 , John P. Wisniewski 25 , Toru 

Yamada 26 , and Motohide Tamura 2,6 


2 


1 College of Science, Ibaraki University, Bunkyo 2-1-1, Mito 310-8512, Japan: 
ttsuka@mx. ibaraki .ac.jp 

2 National Astronomical Observatory Japan(NAOJ), Osawa 2-21-1, Mitaka, Tokyo 181- 
8588, Japan 

^Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, 
Yoshinodai 3-1-1, Sagamihara, Kanagawa 229-8510, Japan 

^Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, 
USA 

5 Department of Physics and Astronomy, The University of Oklahoma, 440 W. Brooks St. 
Norman, OK 73019, USA 

6 School of Science, The University of Tokyo, Hongo 7-3-1, Bunkyo, Tokyo 113-0033, 
Japan 

'Lboratoire Lagrange (UMR 7293), Universite de Nice-Sophia Antipolis, CNRS, Obser- 
vatoire de la Cote d’Azur, 28 avenue Valrose, 06108 Nice Cedex 2, France 

8 Max Planck Institute for Astronomy, Konigstuhl 17, 69117 Heidelberg, Germany 

9 Department of Astrophysical Sciences, Princeton University, Peyton Hall, Ivy Lane, 
Princeton, NJ 08544, USA 

10 Department of Physics and Astronomy, College of Charleston, 58 Coming St., 
Charleston, SC 29424, USA 

11 Subaru Telescope, 650 North A’ohoku Place, Hilo, HI 96720, USA 

12 Universitats-Sternwarte Miinchen, Ludwig-Maximilians-Universitat, Scheinerstr. 1, 
81679 Miinchen, Germany 

13 Exoplanets and Stellar Astrophysics Laboratory, Code 667, Goddard Space Flight Cen- 
ter, Greenbelt, MD 20771 USA 

14 Institute for Astronomy, University of Hawaii, 640 N. A’ohoku Place, Hilo, HI 96720, 


3 


USA 

15 Department of Astrophysics! Sciences, Princeton University, Peyton Hall, Ivy Lane, 
Princeton, NJ 08544, USA 

16 Department of Earth and Planetary Science, The University of Tokyo, 7-3-1 Hongo, 
Bunkyo-ku, Tokyo, 113-0033, Japan 

17 Exoplanets and Stellar Astrophysics Laboratory, Code 667, Goddard Space Flight Cen- 
ter, Greenbelt, MD 20771 USA 

18 Department of Astronomy, Kyoto LIniversity, Kitashirakawa-Oiwake-cho, Sakyo-ku, Ky- 
oto, Kyoto 606-8502, Japan 

19 Hiroshima University, 1-3-2, Kagamiyama, Higashihiroshima, Hiroshima 739-8511, 
Japan 

20 Department of Astrophysics, CAB-CSIC/INTA, 28850 Torrejon de Ardoz, Madrid, Spain 

21 Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, 171-113, 
USA 

22 Institute of Astronomy and Astrophysics, Academia Sinica, P.O. Box 23-141, Taipei 
10617, Taiwan 

23 Astronomical Institute ’’Anton Pannekoek”, University of Amsterdam, Postbus 94249, 
1090 GE, Amsterdam, The Netherlands 

24 Department of Cosmosciences, Hokkaido University, Kita-ku, Sapporo, Hokkaido 060- 
0810, Japan 

25 Department of Astronomy, University of Washington, Box 351580 Seattle, WA 98195, 
USA 

26 Astronomical Institute, Tohoku University, Aoba-ku, Sendai, Miyagi 980-8578, Japan 

2 'The Graduate University for Advanced Studies (SOKENDAI), Shonan International 
Village, Hayama-cho, Miura-gun, Kanagawa 240-0193, Japan 

28 Department of Astronomy and Astrophysics, University of Toronto, 50 St. George Street 


-4- 


Received 


accepted 


Draft version 1 


M5S 3H4, Toronto Ontario, Canada 


- 5 - 

ABSTRACT 

To reveal the structures of a transition disk around a young stellar object 
in Lupus, Sz 91, we have performed aperture synthesis 345 GHz continuum and 
CO (3-2) observations with the Submillimeter Array (~ l /, -3 ,/ resolution), and 
high-resolution imaging of polarized intensity at the AVband by using the Hi- 
CIAO instrument on the Subaru Telescope (Cf/25 resolution). Our observations 
successfully resolved the inner and outer radii of the dust disk to be 65 and 170 
AU, respectively, which indicates that Sz 91 is a transition disk source with one 
of the largest known inner holes. The model fitting analysis of the spectral en- 
ergy distribution reveals an H 2 mass of 2.4 x 10~ 3 M 0 in the cold (T <30 K) 
outer part at 65 < r < 170 AU by assuming a canonical gas-to-dust mass ratio of 
100, although a small amount (> 3 x 10 -9 M 0 ) of hot (T ~180 K) dust possibly 
remains inside the inner hole of the disk. The structure of the hot component 
could be interpreted as either an unresolved self-luminous companion body (not 
directly detected in our observations) or a narrow ring inside the inner hole. Sig- 
nificant CO (3-2) emission with a velocity gradient along the major axis of the 
dust disk is concentrated on the Sz 91 position, suggesting a rotating gas disk 
with a radius of 420 AU. The Sz 91 disk is possibly a rare disk in an evolutionary 
stage immediately after the formation of protoplanets because of the large inner 
hole and the lower disk mass than other transition disks studied thus far. 

Subject headings: stars: pre-main-sequence — (stars:) planetary systems: 
protoplanetary disk — submillimeter — infrared: stars 


1. Introduction 


Nearly all newly-formed stars are surrounded by disks of gas and dust, which provide 
the building blocks of planets (Evans et al. 2009). Thus, studying the structure and 
evolution of these protoplanetary disks provides information about how and when planets 
form. Initially, the disks are optically-thick, producing broadband infrared (IR) emission 
well in excess of the stellar photosphere. The star+disk system has a rather flat IR spectral 
energy distribution (SED) slope (e.g., Lada & Wilking 1984; Lada 1987). By ~ 5 Myr, 
nearly all stars lack evidence of warm circumstellar dust and accretion onto the star (e.g., 
Fedele et al. 2010; Currie & Sicilia- Aguilar 2011). This implies that most protoplanetary 
disks around solar/subsolar-mass stars have disappeared by this time, and the disk material 
has been accreted/dispersed because of processes such as photoevaporation or the disk 
material has been incorporated into planets. 

Transition disks bridge the gap between these endpoints of disk evolution because 
their excess emission at some IR wavelengths is intermediate between those of an optically 
thick disk and a star lacking a disk (e.g., Strom et al. 1989). Although the term “transition 
disk” includes a diverse set of morphologies, a major subset of these objects include those 
with near/mid IR deficits (relative to an optically-thick disk) but optically-thick emission 
at other, typically longer wavelengths, which is indicative of inner holes or gaps (e.g., 
Calvet et al. 2002). IR SED modeling implies that the hole or gap sizes for these objects 
typically range between ~ 1 and ~ 50 AU (e.g., Espaillat et al. 2010; Merin et al. 2010) and 
encloses the planet-forming region in our own solar system. SED modeling has identified 
numerous other transition disks in nearby star forming regions (e.g., Cieza et al. 2010; 
Merin et al. 2010; Currie & Sicilia-Aguilar 2011; Romero et al. 2012; Cieza et al. 2012b). 
Because formation of gas giant planets creates holes or gaps in disks, transition disks with 
SEDs consistent with these features may be excellent laboratories for studying planet 
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formation. 

Although SED modeling can provide indirect evidence of an inner hole or gap, 
high-resolution observations with long baseline (sub-) millimeter interferometers and large 
aperture IR telescopes provide direct evidence of their existence (e.g, Brown et al. 2009; 
Hughes et al. 2009; Andrews et al. 2011; Isella et al. 2012; Cieza et al. 2012a; Mathews et al. 
2012; Thalmann et al. 2010; Hashimoto et al. 2012; Mayama et al. 2012). Since millimeter 
and micron sized dust grains contribute to the majority of millimeter and NIR emissions, it 
is important to compare the detailed spatial distributions of these tracers for understanding 
the physics of the transition disks. 

Sz 91 is an M0. 5 young star surrounded by a transition disk located in the Lupus 111 
molecular cloud ( d = 200 pc; Comeron 2008). The stellar position, (aj 2 ooo, <5j20oo)? is (16 h 
07 m 11?6, —39° 03' 47'.'2). The stellar mass and age have been estimated to be 0.49 M 0 and 
5 Myr, respectively (Hughes et al. 1994). The Ha line width at 10 % of the peak has been 
measured to be 283 km s -1 , suggesting accretion with 10” 10 M 0 yr" 1 (Romero et al. 2012). 
Sz 91 is classified as a wide binary with a separation of 9"(Melo 2003), corresponding to 
1800 AU. However, such a wide binary system is unlikely to disturb the disk evolution of 
a host star (Kraus et al. 2012; Harris et al. 2012). Moreover the difference between their 
proper motions implies that these sources are not co-moving (Roeser et al. 2010). We 
therefore treat the Sz 91 disk as a circumstellar disk around a single star system throughout 
this paper. 

The SED of Sz 91 is characteristic of the transition disk. The following remarkable 
features of this source can be identified in the SED: no significant IR excess with a spectral 
index from K s to 24 //m of —2 (i.e. , class III in the IR categorization), presence of a large 
dip of approximately 20 //in, and very steep flux density rising between 24 and 70 /mi 
(Evans et al. 2009; Romero et al. 2012). Romero et al. (2012) categorized Sz 91 as a giant 


planet-forming disk on the basis of the following features: very steep increase in flux density 
at 24 /iin, which indicates a sharp edge at the inner radius of the disk; a clear sign of mass 
accretion onto the central star; a relatively massive disk (~5x 10- 3 M 0 ). 

Although most transition disks with inner holes have flux deficits restricted to 1-10 
/irn and optically thick emission at longer wavelengths, Sz 91 differs from them because 
its flux deficit extends to significantly longer wavelengths (24 /irn) , and it has strong 
far-lR to millimeter emission. The total flux density of the 870 /irn (345 GHz) continuum 
emission of Sz 91 has been estimated to be 34.5±2.9 mJy (Romero et ah 2012). Although 
the 1.3 mm continuum emission had not been detected previously (Nuernberger et al. 
1997), the recent wide held imaging survey with the AzTEC receiver on the Atacama 
Submillimeter Telescope Experiment (ASTE) has clearly detected a 1.1 mm flux density 
of 27.2±6.0 mJy (Kawabe et al. 2014). These submillimeter flux densities imply that the 
Sz 91 ’s disk is substantially massive with respect to those of other class 111 sources in 
nearby star forming regions (Andrews & Williams 2005, 2007); the flux densities roughly 
correspond to 2-5 xlO" 3 M 0 if the canonical gas-to-dust mass ratio of 100 is assumed. Such 
an object having significant submillimeter flux density, but no NIR excess appears to be 
at the transition phase from class II to III and is quite rare (~5 % Andrews & Williams 
2005, 2007). Although a recent high- resolution imaging survey for transition disk objects 
(Andrews et al. 2011) has successfully revealed a large inner hole in the disk up to a radius 
of ~ 70 AU, there is selection bias for the most massive transition disks. Because the 
submillimeter flux density of Sz 91 is the lowest among the transition disk objects studied 
thus far, considering the distance to the source (Fig. 10 in Andrews et al. 2011), Sz 91 is 
a crucial target for investigating the disk evolution from the perspective of a variety of 
submillimeter flux densities, i.e., the disk mass. However, the disk structure has not been 
resolved thus far, and high-resolution imaging observations are urgently required. 
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In this paper, we present the first sub-arcsecond resolution images toward Sz 91 at 
submillimeter and near infrared wavelengths. In section 2, observational parameters and 
calibration for the obtained data are described. In section 3, high- resolution images are 
shown. In section 4, the detailed disk structure is discussed on the basis of model fitting to 
the SED and the CO (3-2) line profile. The unique nature of Sz 91 is determined from a 
comparison with other transition disks. Finally in section 5, the key results of this study 
are summarized. 


2. Observations and Data Reduction 

2.1. High-resolution Imaging of the 345 GHz Continuum and CO(3— 2) 

Emission 

Interferometric observations of the 345 GHz continuum and CO (3-2) emission toward 
Sz 91 were conducted in 2010 with the Submillimeter Array (SMA; Ho et al. 2004), which 
comprises eight 6 m antennas located atop Mauna Kea in Hawaii. The observations 
were performed in the compact configuration with baselines ranging from 6 to 70 m. For 
the continuum emission, higher resolution data were also obtained in 2010 in the very 
extended configuration with baselines up to ~ 500 m. The field of view of SMA was 
approximately 30" in full width at half maximum (FWHM). The emission was detected 
using double side-band superconductor-insulator-superconductor receivers with a local 
oscillator frequency of i'lo = 340.755 GHz. Both the upper and lower sidebands data were 
used for the continuum, resulting in an 8 GHz bandwidth in total. The channel spacing of 
the CO(3-2) line was set to 0.8125 MHz. 

The amplitude and phase of the array system were calibrated by observations of 
quasars, J1604— 446, J1626— 298, and J1517— 243, in a cycle with 12.5 min on a target and 
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6 or 8 min on a quasar. The response across the observed passbands was determined by 60 
min observations of quasar, 3C 273. Absolute flux calibration was achieved by observing 
Titan at the beginning and end of each night. 

We created the continuum images by combining the very extended and compact 
configuration data, whereas the CO maps were created from the compact configuration data 
only. The UV data were edited and calibrated using MIR, an IDL-based software package. 
We used Astronomical Image Processing System (AIPS) for imaging procedures, including 
deconvolution by the CLEAN algorithm and restoration with a synthesized beam. Both 
continuum and line emission maps were created with natural weighting in the visibility 
plane, producing synthesized beam FWHMs of 1"A x l'.'O at a position angle (PA) of — 9?8 
and 4'.'6 x 1"8 at a PA of — 4?0, respectively. 


2.2. High-resolution Polarized Intensity Imaging at K s band with the Subaru 

Telescope 

iC-band (2.15 //m) linear polarized intensity (PI) images of Sz 91 were obtained 
with the high-contrast imaging instrument (HiCIAO; Tamura et ah 2006) combined with 
dual-beam polarimetry from the 8.2 m Subaru Telescope in May 2012. The observations 
were conducted under the program SEEDS (Strategic Explorations of Exoplanets and Disks 
with Subaru; Tamura 2009). The adaptive optics system (A0188: Hayano et al. 2004) 
provided a limited diffraction and mostly stable stellar point spread function (PSF) with 
a FWHM of 0"06 in the K s band. Polarization differential imaging (PDI) is a powerful 
technique used to reveal the structure of a dusty disk in very close proximity to a star 
(e.g., Thalmann et al. 2010; Hashimoto et al. 2011, 2012; Tanii et al. 2012; Muto et al. 
2012; Kusakabe et al. 2012; Mayama et al. 2012). We employed the PDI mode, combined 
with the angular differential imaging mode, in which the field of view and the pixel scale 
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were 10" x 20" and 9.5 mas pixel -1 , respectively. Half-wave plates were placed at the 
four angular positions of 0°, 45°, 22?5, and 67?5 in sequence with a 30 sec exposure per 
wave plate position. Image Reduction and Analysis Facility (IRAF 1 ) software was used 
for all data reduction as following the methods of Hashimoto et al. (2011) and Tanii et al. 
(2012), and the Stokes Q and U parameter images were created. We calculated the PI as 
\J Q 2 + U 2 and the polarization vector angle as 0.5 x arctan (U/Q) with the 3.4.0 version of 
the Common Astronomy Software Applications package. The final PI image was created 
from the smoothed Q and U images by a Gaussian function; thus, the effective resolution of 
the image became ~ 0"25. 


3. Results 

3.1. 345 GHz Continuum Emission Map 

We created a 345 GHz continuum emission map of Sz 91 with all UV data as shown 
in Figure 1(a). The map shows a strong continuum emission peak at the stellar position. 
The total flux density is measured to be 32.1±3.6 mJy, which agrees well with the flux 
density of 34.5±2.9 mJy, which is the previous 345 GHz measurement with the 12 nr single 
dish APEX telescope obtained by Romero et al. (2012). The 345 GHz emission is resolved 
by the F/4 x l'.'O beam at a PA of — 9?8. The beam- deconvolved size of the emission is 
measured to be (l'/7±0'. / l) x (O'.TiO^) with a PA of 169?0±3?0 from a 2-D Gaussian fitting, 
corresponding to (340±20) x (140±40) AU at a distance of 200 pc. This size is comparable 
to that of a typical protoplanetary disk. The 345 GHz emission is probably thermal dust 

1 IRAF is distributed by National Optical Astronomy Observatory, which is operated by 

the Association of Universities for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 
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emission from the disk around Sz 91 because the flux densities from 350 to 1100 pm 
monotonically decrease with a spectral index of ~ 2 (Figure 5), which is comparable to 
other T Tauri stars (Andrews & Williams 2005, 2007). 

Moreover, we created a higher-resolution 345 GHz image using the UV data over the 
UV range from 50 to 600 kA as shown in Figure 1(b). The emission peak at the 4.5a level 
can be found to the north of the star. The peak position is measured to be (aj 2 ooo =16 h 
07 m 11!6, <5j2ooo = —39° 03' 4778) and is shifted from the stellar position by 0743 toward 
the direction at a PA of 1?4. The lack of a peak at the stellar position implies that the 
inner part of the dust disk is depleted or cleared as expected from the SED. Therefore, the 
continuum emission in Figure 1(b) probably originated from the innermost part of the dust 
disk. Assuming that the peak position in figure 1(b) represents the inner edge of the dust 
disk, the inner radius of the Sz 91 transition disk is estimated to be 86±25 AU. The total 
flux density of this compact component is measured to be 13.2±3.1 mJy, which is 40 % of 
the total flux density in Figure 1(a). 


3.2. CO(3— 2) Molecular Line Emission Maps 

Figure 2 shows velocity channel maps of the CO (3-2) line. The significant CO (3-2) 
emission is detected in the velocity range from 1.9 to 7.2 km s -1 in Vlsr- It is clear that 
the emission at Vlsr =1.9-5. 8 km s -1 is concentrated at the stellar position. In addition, 
a velocity gradient appears to be present along the north-south direction in this velocity 
range: the emission is located mainly on the north side of the star at Vlsr =2. 6-2. 9 km 
s~ x and at the south from 4.0 to 5.1 km s -1 . Notably, the central velocity range from 3.3 
to 4.0 km s” 1 includes the radial velocity of Sz 91 (Vlsr =3.87 km s -1 , converted from 
Vhoiio = —1.57 km s _1 ; Melo 2003). In the range of 4. 4-5. 4 km s _1 , we detected spatially 
extended emission near the star (4.4 and 4.7 km s" 1 ) and at the south side of the held of 
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Fig. 1. — CLEANed maps of the 345 GHz continuum emission created by all UV data (a) 
and 50-600 kA (b). The contours begin at ±2cr with the intervals of 2a, where la =2.1 and 
2.8 mJy beam” 1 in (a) and (b), respectively. The red cross indicates the stellar position in 
each panel. The synthesized beam sizes in half power beamwidth (HPBW) are shown at the 
bottom- left corner in each panel: 1 "A x l'.'O with a position angle (PA)=— 9?8 and 0"7 x 0'.'3 
with a PA=13?0 for (a) and (b), respectively. The blue box in the left panel indicates the 
area of the right panel. 
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view (4. 7-5.4 km s' 1 ). The origin of the extended components is most probably an ambient 
cloud because the systemic velocity and the velocity width of the main cloud condensation 
of Lupus 111 have been measured to be 4.1 and 1.2 km s' 1 , respectively (Hara et al. 1999). 
We also detected the extended emission in the 6. 1-7.2 km s' 1 channels. These components 
are shifted by > 3" toward the south-west direction from the star and also possibly originate 
from the spatially extended ambient cloud. Hereafter, we focus on compact components 
toward the star that should have originated from the gas disk; we do not discuss the 
extended components in this paper. 

To clearly see the gas disk emission and its velocity gradient, we created total integrated 
intensity (1. 5-6.1 km s' 1 ) and first moment maps of the CO(3-2) line, as shown in Figure 
3. From these maps, we successfully detected the centrally concentrated emission and the 
velocity gradient roughly along the north-south direction, suggestive of the rotating gas 
disk. The total integrated intensity of CO(3-2) in the area above the 3cr noise level is 8.55 
Jy km s' 1 . The emission peak is shifted from the stellar position by ~ 1" toward the north 
direction, possibly because the red-shifted emission is partially resolved-out owing to the 
contamination by the ambient cloud. 


3.3. Polarized Intensity Image at K s band 

The high-resolution PI and polarization vector angle maps at the AT s -band are shown 
in Figures 4(a) and (b), respectively. The presented images are smoothed and the effective 
resolution of the images is ~ O'.^S. A crescent-like emission region was detected around 
Sz 91, which is elongated from the south to the north via the west side of the star. The 
polarization angles are nearly perpendicular to the radial directions from the central star, 
indicating that the J\ s -band emission probably originated from the scattered light at the 
inner part of the dust disk. The crescent-like emission suggests that this part is the near 


15 



Fig. 2. — CO(3-2) velocity channel maps with 0.35 km s -1 resolution obtained by SMA. 
The central LSR velocity of each panel is shown at the top in km s -1 . The contours start 
at ±2cr with intervals of 2cr, where lcr = 289 mJy beam -1 . The ellipse at the bottom left in 
the top left panel shows the synthesized beam size of 4"6 x l'/8 with a PA of — 4?0 in HPBW. 
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Fig. 3. — Total integrated intensity map of the CO(3-2) emission over a Vlsr range of 1. 5-6.1 
km s -1 (contours) superimposed on the intensity- weighted first moment map (color). The 
contours start at ±3cr with intervals of 3cr, where la =450 mJy km s -1 beam -1 . The color 
bar on the right-hand side of the panel shows the centroid velocity in m s -1 . The red cross 
indicates the stellar position. The synthesized beam size in HPBW is shown at the bottom 
left corner: 4'.'6 x 1''8 with PA= — 4?0. 
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side of the disk if we assume that forward scattering is dominant, as is the case in Mie 
scattering. Substantial emission dips appear at the north and south sides of the software 
mask, possibly due to the inner hole structure of the dust disk. The PSF subtraction 
process is a primary factor to cause an artificial systematic error near the star, and the 
error becomes a systematic emission via the square root operation of stokes Q and U. The 
symmetric distribution of the dips indicates the existence of an inner hole in the dust disk. 

We also detected a bright region near the western edge of the software mask, i.e. , a 
photometrically reliable emission. However, the origin of this emission remains unclear. 
Because the PSF shows a substantially symmetric distribution, it is unlikely that the bright 
emission would appear at the only one-side. Because it extends to the edge of the software 
mask, higher resolution and higher contrast observations close to the star are required to 
determine the origin of the bright emission. 

The existence of the inner hole structure expected from the PI image strongly supports 
the results of the submm images. The higher-resolution 345 GHz image is overplotted in 
Figure 4(a). It is clear that the PI is distributed at the innermost part of the submm 
emission. For example, along PA=0°, the 345 GHz emission peaks at C//45 from the star 
whereas the PI peaks at 0 , /32. These facts suggest the existence of the hole structure at the 
inner part of the disk and that the PI originates from the innermost part of the transition 
disk, i.e., the inner edge. 

To quantify the shape of the inner edge of the dust disk, we fit an ellipse to the 
crescent-like emission. Table 1 lists the parameters of the best-fit ellipse, and the fitted 
ellipse is shown in Figure 4(a) by the yellow dotted line. The center position of the ellipse 
coincides with the stellar position. The PA of the ellipse, i.e., the direction of the disk 
major axis, is along at 17?5 ± 17?7, which is consistent with the PA of the emission dips 
in the crescent-like emission. The PI at the disk minor axis (PA~ 120°) appears to be 


lower than that at the major axis. This result is probably due to lower polarization degree 
along the minor axis where the scattering angle at the disk surface is deviated from 90° 
(McCabe et al. 2002). 


4. Discussion 

4.1. Spatially Resolved Disk Structure 

Both the high- resolution subnnn and NIR images clearly resolved the disk structure 
around Sz 91. In this section, we estimate the inner radius, outer radius, and inclination 
angle of the disk directly from our images. These parameters are adopted in the model 
calculations in sections 4.2 and 4.4. 

The inner radius of the dust disk can be estimated from the large inner hole revealed 
in the AVband PI image. Assuming that the AVband emission is the scattered light at 
the inner edge of the disk, the best-fit ellipse in the AVband image shows the inner radius 
to be 65±4 AU. There are two possibilities to explain the scattered light. One is forward 
scattering at the surface of the inner part of the disk, and the other is scattering at the 
inner wall of the disk such as that in the case of LkCa 15 (Thalmann et al. 2010). We 
can not judge which case is more likely to explain the NIR emission from the Sz 91 disk. 
However, in either case, it probably originates from the innermost part of the disk, and 
our estimation of the inner radius seems to be reasonable. The high-resolution 345 GHz 
image also provides a rough estimate of the disk inner radius. Although the 345 GHz 
image in Figure 1(b) does not clearly show the presence of the inner hole, it is clear that 
the 345 GHz flux density becomes weak toward the central star, which is suggestive of dust 
depletion near the star. The angular distance between the 345 GHz emission peak and the 
central star can be interpreted as the disk inner radius of 0"43 (= 86 AU), which is roughly 
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Fig. 4. — (a) Polarized intensity (PI) image obtained with the Subaru telescope in the 
K s band (color) where the 345 GHz continuum emission map (contour) in Figure 1(b) is 
superposed. The stellar position is shown by the red cross. The best-fit ellipse is indicated 
by the yellow dashed line. The inner area of 0"25 diameter (filled circle) is photometrically 
unreliable owing to the point spread function (PSF) subtraction process and is masked, (b) 
Polarization vector angle map superimposed on the PI image in the left panel. The vector 
directions indicate angles of polarization. The vector’s lengths are arbitrary. 
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consistent with the above estimate of 65 AU. In the following discussion, we adopt the disk 
radius R in of 65 AU. 

The outer radius of the dust disk can be estimated from the beam-deconvolved size 
of the full-UV 345 GHz map shown in Figure 1(a). Assuming a geometrically thin disk, 
the beam- deconvolved size of (l"7±0 , .'l)x(0"7±0"2) corresponds to the outer radius R out of 
170±20 AU at the distance of 200 pc. The gas disk is also resolved in CO (3-2) (Figure 3). 
The beam-deconvolved size is measured to be (3'. , l±0 , . / 3)x(l'. , 6±0 , . , 7) at a PA of 15±ll°from 
a 2-D Gaussian fitting, corresponding to R out =310±60 AU. However, the red-shifted 
gas is probably contaminated by the ambient cloud gas, and the center position of the 
beam-deconvolved disk is shifted by 0 % from the star to the north. Thus we consider 
the largest distance of 2 , . , 1 from the stellar position as the outer radius of the gas disk: 
R out =420 AU. Because the actual intensity distribution of each emission is expected to have 
a power-law form, these outer radii derived from the Gaussian fitting should be regarded as 
lower limits. 

The outer radius of the gas disk is much larger than that of the dust disk. The 
discrepancy in radial extension between the dust and gas emissions is frequently observed 
for other protoplanetary disks, which can be explained by an exponential decrement in 
the surface density (e.g., Hughes et al. 2008; Panic et al. 2009; Isella et ah 2010) or by the 
radial drift of large grains (Andrews et ah 2012). In addition, the possible north-south 
asymmetry in the dust continuum emission, suggested in Figure 1(b), could have resulted 
from the azimuthal drift of large grains (Birnstiel et ah 2013; van der Marel et ah 2013) 
associated perhaps with a perturbing body orbiting within the central hole of this disk. 

The disk inclination angle can be derived from the ratio of the major to minor 
axis lengths of the best-fit ellipse in Figure 4(a): 40±15° from face-on. Moreover, the 
beam-deconvolved size of the 345 GHz continuum emission in Figure 1(a) provides an 
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inclination angle of 66±5°. The latter value of 66°, however, is possibly influenced by the 
disk flaring because the distribution of the 345 GHz emission is biased to the outer part 
of the disk. Thus, we adopt the AVband value of 40° as the disk inclination angle in this 
study. In fact, for the CO(3-2) spectrum, the model calculation in section 4 prefers the 
smaller value of 40°. 


4.2. Disk Parameters Derived from SED fitting 

To deduce other parameters of the dust disk, we performed least-square fitting to the 
SED on the basis of the power-law disk model. Figure 5 shows the SED of Sz 91 that 
includes our 345 GHz flux density and a recent data set by WISE. As Romero et al. (2012) 
reported, the SED is characterized by large dip around 20 pm, sharp rise from 20 to 70 pm, 
and significant (sub-)millimeter emission. 

We introduced the two components of a cold disk and a small amount of hot dust 
inside the disk as follows. First, we applied the model for the cold disk, which comprised 
the usual power-law disk and blackbody stellar emissions, to all of the SED data except 
that at 10-30 pm, because the model could not reproduce all the data points including the 
three points at 10-30 pm. The power-law disk model has surface density and temperature 
radial distributions of E(r) and T(r), respectively, in a power- law form (Kitamura et al. 
2002; Tsukagoshi et al. 2011): 


E(r) 

-*(£) 

(i) 

T(r) 

/ r \ 

(2) 

= T -Ui) • 


where r is the radial distance, Ei n and T m are the surface density and temperature at the 
inner radius, 7?, in , and p and q are power-law indexes. The gas-to-dust ratio was assumed to 
be 100 and the extinction of the stellar light because of the interstellar dust was corrected 
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by the Ay value toward the star and the dust mass absorption coefficient (Figure 1 of 
Adams et al. 1988). Table 2 lists the fixed parameters in the SED fitting. The stellar 
mass, the effective temperature of the star, and the visual extinction toward the star are 
from Hughes et al. (1994). The inner and outer radii and inclination angle of the disk 
were determined in section 4.1. We here adopted the values of 1.5 and 0.5 as p and q, 
respectively, which are the same as those of the minimum mass solar nebular (Hayashi 
1981). These values were selected because the SED fitting is known to be insensitive to the 
power-law index of p and it is difficult to resolve the parameters for the temperature profile 
(T and q) from only the longer wavelength data (> 30 /mi). Notably, the lower limit of 
the disk temperature was set to be 10 K, which is the typical temperature of the Lupus III 
cloud (Vilas-Boas et al. 2000). The stellar radius, A*, T in , E in , and the power- law index of 
the dust mass opacity coefficient j3 (k v = 0.1 x (y / 10 12 Hz)^ cm 2 g -1 : Beckwith et al. 1990) 
were treated as the free parameters in the SED fitting. 


After the power-law disk model fitting, we introduced the additional hot component 
inside the disk to reproduce the observed flux densities at ~ 20 /mi. The presence of 
the inner hot component is supported by the sign of mass accretion derived from the Ha 
emission line (Romero et al. 2012). Because no information was available on the structure 
of this component, we simply assumed a gray body with a temperature of T c , a column 
density of £ c , and a solid angle of fl c . Its flux density S\ at wavelength A is written by 


S\ = 


2 he 2 1 

ex P(Aife) - 1 


(l_ e — K „S C ) x Qc ^ 


( 3 ) 


where h is the Plank constant, c is the speed of light, and k-g is the Boltzmann constant. We 
applied the same dust mass opacity coefficient, k u , as that in the cold disk. The parameters 
of T c , E c , and fl c were treated as free. In the fitting, the upper limit of fl c was set to be 
~ 6.0 x 10' 12 str, corresponding to the solid angle of the inner hole seen in the NIR image. 


The best-fit model SED of the cold disk reproduced the longer wavelength data 
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(> 30 //in) effectively, as shown in Figure 5. The best- fit parameters of the cold disk are 
summarized in Table 3. The (3 value of 0.5±0.1 is significantly smaller than that in the 
diffuse interstellar medium of ~ 2 (Draine & Lee 1984), suggesting that the dust growth 
occurs in the disk (Miyake & Nakagawa 1993). By adopting this (3 and introduced by 
Beckwith et ah (1990), who assumed a 100:1 mass ratio between gas and dust, the disk 
mass is derived to be (2.4±0.8)xl0~ 3 M & . This value is significantly higher than those of 
other class III sources in nearby star forming regions (Andrews & Williams 2005, 2007); 
however, the value is lower than those of most of the transition disks studied thus far 
(Andrews et al. 2011), even if we consider the difference in by a factor of ~1.7 at v = 340 
GHz between this study and the previous studies in which (3 is set to be 1. 

To verify the fitting result, we also created a continuum image for the best-fit model 
and compared it with the observations, as shown in Figure 6. The model image essentially 
agrees with the observations, but there remains a difference of at most 4cr; the negative 
residual extends in the east-west direction, which may be due to the asymmetry of the disk. 

We next attempted to reproduce the observed flux densities at 10-30 // m by adding 
the contribution from the hot component. The best-fit parameters of T c , E c , and to 
reproduce all of the SED data were searched in the reduced y 2 maps by manually changing 
the initial values of T c , E c , and with 6000 runs. 

From the calculations, we determined that there are two distinct regions in the E c -ff c 
plane where the reduced y 2 takes its local minimum values, which indicates a strong 
coupling between E c and fi c . One is the region in which E c < 5 x 10 5 g cm 2 and 
b2 c > 4 x 10” 13 str, and the emission of the hot component is substantially optically thin. 
In this region, the best-fit T c typically converges to be 186 K and the best-fit E c is inversely 
proportional to f2 c . The product of E c x is a constant value, which provides the best-fit 
mass of 3 x 10~ 9 M 0 for the hot component. The other is the region in which E c > 30 
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g cm -2 and the emission is substantially opticall thick. In this region, the best-fit T c is 
typically 172 K and X c is independent of fl c . For the hot component, we can obtain the 
best-fit H c = lx 10“ 16 str and the lower limit mass of 6 x 10“' M 0 . Although the reduced 
X 2 values are slightly lower in the former case (Ay 2 = 0.3), the difference between the two 
regions is not significant. We therefore conclude that there are two possible origins of the 
hot component: the optically thin gray body emission and the optically thick black-body 
emission. However, the optically thin condition is the unlikely case because the mass of the 
hot component is too small given the mass accretion rate of 10“ 10 M 0 yr _1 ; all the hot 
component would disappear in only 30 yrs. 

Notably, the contribution of the hot component is mainly restricted by only three data 
points at MIR, from 12 to 25 /irn: the upper limit of the temperature is determined by the 
slight excess emission at 12 /irn and the lower limit is limited by the decrement between 22 
and 24 /irn. A refined model will be required when more data at MIR and FIR are obtained 
by further observations. 


4.3. Origin of the Hot Component in the Inner Hole of the Disk 

The SED model analysis indicates that the presence of a hot component inside the 
dust disk with a temperature of ~180 K, which is responsible for the SED peak at ~20 /irn. 
There are two possible origins for such a hot component: a localized self-luminous emitting 
body (i.e. , circumplanetary disk) or an inner warm structure of the disk (Wolf & D’Angelo 
2005). In this section, we estimate the size and mass of the hot component for both cases. 

In the case of the circumplanetary disk, we can estimate the radius from f2 c because 
it is usually optically thick (D’Angelo et al. 2003). Assuming the circumplanetary disk is 
parallel to the parent circumstellar disk, the solid angle of 1 x 10“ 16 str corresponds to the 
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Tabic 1. Parameters of the best-fit ellipse in Figure 4 


Parameter 

Best-fit value 

Error 

AR.A. 1 (") 

0.03 

0.04 

ADec. 1 (") 

0.03 

0.02 

Major axis (") 

0.33 

0.02 

Minor axis (") 

0.25 

0.08 

Position angle (°) 

17.5 

17.7 


1 Offset from the stellar position. 


Table 2. Fixed parameters in the SED fitting 


Parameter 

Fixed Value 

M* (Mq) 

0.49 1 

Ay (mag) 

2.0 1 

Rin (AU) 

65 

Rout (AU) 

170 

i (degree) 

40 

P 

1.5 

q 

0.5 


1 Hughes et al. (1994) 
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1 10 100 1000 

X (pm) 

Fig. 5. — Spectral energy distribution of Sz 91. The filled circles indicate observed 
flux densities. All data except the 345 GHz flux density were compiled from previ- 
ous studies, including the NOMAD catalog (Zacharias et al. 2005), the 2MASS point 
source catalog (Cutri et al. 2003), the Spitzer IRAC and MIPS photometry (Evans et al. 
2009), the WISE all-sky data release (Cutri & et al. 2012), the AKARI FIS all-sky sur- 
vey point source catalog (Yamamura et al. 2010), the IRAS point source reject cat- 
alog (Infrared Astronomical Satellite IRAS), the AzTEC/ASTE 1100 /irn photometry 
(Kawabe et al. 2014), and the Herschel PACS/SPIRE FIR flux densities which we derived 
from the HSA Science archive data. The best-fit SED and the contribution of the hot com- 
ponent in the optically thin case are shown by the black and red lines, respectively, and by 
the gray and orange lines in the optically thick case. The blue dashed line indicates the 
stellar contribution. 
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Table 3. Best-fit parameters in the SED fitting 


Parameter 

Value 

Stellar parameters 

R* (Rq) 

1.29±0.03 

T* (K) 

4148T55 

L* (Rq) 

0.49T0.02 

Cold outer disk 

Tin (K) 

32.5T3.9 

Sin 1 * (g cm -2 ) 

0.67T0.03 

p 

0.54=0.1 

Mdisk 1 (10“ 3 M 0 ) 

2.44=0.8 


1 The surface density and the 

disk mass are shown in the 

gas+dust density and mass by as- 
suming the gas-to-dust mass ra- 
tio of 100. 


radius of 0.3 AU, or 64 R Q . On the other hand, the Hill radius of a putative planet around 


Sz 91 is expressed by 


= 1.1 x 


10 AU VI Mj 


m 


1/3 


[AU], 


( 4 ) 


where m and a are the mass and orbital radius of the planet, respectively. The radius of the 
circumplanetary disk inferred from the SED fitting is significantly smaller than -Rhm at 3-65 
AU from the star, which is consistent with the theoretical expectation for a circumplanetary 
disk (e.g., Tanigawa et al. 2012). However, our ADI observation with the Subaru telescope 
could not examine the presence of such a companion planet due to a low rotation angle of 


~ 16°. 


For the second possibility, because the hot component could be fitted with a single 
temperature, it must be confined to a narrow width in radius. If we extrapolate the 
temperature distribution of the outer disk determined in the SED fitting, the temperature 
of the hot component in the optically thick case (172 K) corresponds to a radius of 2.3 
AU. The best-fit fi c = 1 x 10 -16 str corresponds to a ring width of 0.01 AU, which is 
significantly narrow with respect to the ring radius. In contrast, the optically thin condition 
is unlikely because the solid angle is > 4 x 10 -13 str which corresponds to a ring width 
of >53 AU. Therefore we conclude that the other implication of the hot component is 
the optically thick ring at 2.3 AU, whose total mass is at least 6 x 10 -7 M 0 . Such an 
example of the optically thick ring around the transition disk has also been reported in 
RX J1633. 9-2442 (Cieza et al. 2012a). The inner structure of Sz 91 may be similar to that 
of RX J1633. 9-2442, whereas Sz 91 exhibits a larger inner hole and a lower disk mass. 
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4.4. Velocity Structure of the Gas Disk: Model Calculation and Comparison 


The CO(3-2) image shown in Figure 3 suggests the presence of a rotating gas disk 
around the star. To reveal the disk rotation in detail, we calculated model spectra of 
CO (3-2) with a simple power-law disk model according to that reported by Kitamura et al. 
(1993), and we compared the results with the observed CO (3-2) profile. 

The following parameters of the model disk were estimated from the observed images 
and the SED fitting: the inner and outer disk radii of 65 and 420 AU, respectively; the 
temperature distribution of T(r) = 32.5 x (r/i? in ) -0 ' 5 K, and the surface density distribution 
of E(r) = 0.67 x (r/i7 in )" 1,5 g cm -2 . The hydrostatic equilibrium is assumed along the 
vertical direction and the density distribution, p(r,z), is therefore expressed by 


Here, G is the gravitational constant, M* is the stellar mass, p is mean molecular weight, 
and me is the mass of H atom. The density at the midplane is expressed by 



The gas motion is assumed to be Kepler rotation whose velocity field is written by 


with CO (3— 2) Profile 



( 5 ) 


where H[r) is the scale height given by 



( 6 ) 



( 8 ) 


The fractional abundance of CO with respect to H 2 is assumed to be 9 x 10“ 5 , which 
corresponds to a typical interstellar value (e.g., Schloerb & Snell 1984; Irvine et al. 1985). 
Because the density of the disk is at least 10 6 cm -3 at the midplane which is significantly 
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higher than the critical density of CO(3-2) (~ 10 4 cm -3 ), we assume local thermodynamic 
equilibrium (LTE). Notably, the hot component in the dust hole was not included because 
its contribution was negligible (<lx 10 -4 Jy). 

Figure 7 shows the calculated CO spectra superimposed on the observed CO spectrum 
integrated over a 5" x 5" box centered at the stellar position. We noted that the peak 
intensity of the calculated CO spectrum differed from the observed value. However, the 
discrepancy is not significant because the disk emission around the stellar LSR velocity 
is probably contaminated by the ambient cloud emission. Therefore, we focused on the 
emission at the blue-shifted side (Vlsr < 2.9 km s _1 ). The model profile in the i = 40° case 
agrees with the observed data, confirming the validity of the disk parameters derived from 
the dust disk. Notably, the inclination angle of 66° determined by the 345 GHz image did 
not fit well the CO line shape, indicating that the inclination of the dust disk estimated 
from the NIR scattered light is more plausible. 

The presence of the gas inside the inner edge of the dust disk could not be confirmed 
from our data set. Although we calculated the model spectrum for R m =0 AU as shown in 
Figure 7, the difference is within the la uncertainty. The presence of the gas disk in the 
inner hole is supported by the fact that Sz 91 shows the mass accretion onto the star, and 
thus, higher spatial resolution imaging is required to reveal the inner structure of the gas 
disk. 


5. Summary 

We present the results of the aperture synthesis 345 GHz continuum and CO (3-2) 
line emission observations with SMA toward a transition disk object in Lupus, Sz 91. 
Furthermore, a high resolution image of the polarized intensity at the K s band obtained 
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with the Subaru Telescope is also presented. The transition disk around Sz 91 has been 
directly resolved and imaged in this study. The disk parameters are derived and the 
structure and evolutional phase of the Sz 91 transition disk are discussed. The main results 
of our observations are summarized in the following points: 

• Our high-resolution imaging revealed a dust disk around Sz 91 with inner and outer 
radii of 65 and 170 AU, respectively, and an inclination angle of 40°. Furthermore, 
the Kepler rotating gas disk with a radius of 420 AU was imaged in the CO (3-2) line. 

• Model analysis of the SED of Sz 91 was performed by using a simple power-law 
disk model. We determined that the observed SED can be reproduced well by the 
combination of a cold disk and a hot component in the inner hole of the disk. The 
total H 2 mass of the cold disk is estimated to be 2.4 x 10 -3 M 0 if the canonical 
gas-to-dust mass ratio of 100 is adopted. The disk mass is significantly higher than 
those of other class III sources in nearby star forming regions; however, the disk mass 
is one of the lowest masses among the currently known transition disks. 

• We determined that the hot component can be expressed by a single temperature 
gray body of ~ 180 K. Although the hot component could not be resolved by our 
observations, its origin is either a localized self-luminous emitting body (i.e. , a Jovian 
mass protoplanet with a circumplanetary disk) or an optically thick ring in the inner 
hole of the disk at 2.3 AU. 

• Our results confirm the previous results such that the disk structure of Sz 91 is 
consistent with that of an ongoing giant planet forming disk. In particular, the 
relatively large inner hole and lower disk mass indicate that the transition disk of 
Sz 91 is probably in a stage of nearly completing planet formation. Sz 91 will be a 
crucial target for investigating the evolution of transition disks and the planetary 
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formation process. In the near future, our proposed study with the Atacama Large 
Millimeter/submillimeter Array (ALMA) will provide a new insight into the planet 
formation process. 
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supported by JSPS KAKENHI Grant Numbers 24103504 (T.T.) and 23103004 (M.M.). 
J.C. gratefully acknowledges support from NSF grant AST-1009203. The Submillimeter 
Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 
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Fig. 6. — Comparison of the observed data with the best-fit model. The 345 GHz continuum 
image, model image, and residual image (data— model) are shown from left to right. The 
contour starts at ±2cr, where la =2.1 mJy beam -1 . The solid and dashed lines indicate 
positive and negative, respectively. 
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Fig. 7. — CO (3-2) spectrum integrated over a 5" x 5" area centered at the stellar position 
(black line). The red line shows the calculated profile by a disk model with an inner radius 
of 65 AU and an inclination angle of 40° from face-on (see text for the detailed description 
of the model). The blue line indicates the case of the inner radius of 0 AU. The green line 
indicates the model profile of the disk with a 65 AU inner radius and an inclination of 66°. 
The CO emission in the velocity range from 4 to 7 km s _1 in Vlsr is probably contaminated 
by the ambient cloud emission, as shown by the gray line at the top of the panel. 
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